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Abstract. The Hubble Deep Field (North) and the flanking fields are used investigate the occur- 
rence of multiple weak lensing deflections along the line of sight in relatively deep imaging data 

■ (zicns ~ 0.6, Zsourco ~ 1-2). Ray tracing simulations of galaxy-galaxy lensing in the HDF-North 
show that proper inclusion of multiple weak deflections is important for a correct prediction 
of the net shear for most sources, and for a given source redshift the number of multiple weak 
deflections is largely insensitive to the cosmography. The effects of multiple weak deflections on 
the magnitude of the weak lensing signal are, of course, strong functions of the adopted halo 

^ ■ parameters. Independent of the halo parameters, however, the closest lens to a source (in pro- 

■ jection on the sky) is not the strongest lens in the case of more than 50% of the sources which 
acquire a net shear of 7 0.01. In addition, multiple weak deflections result in a tangential 
shear about the lens centers that is greater than the tangential shear that would occur if source 
galaxies were lensed solely by the closest lens. Further, multiple weak deflections give rise to 

I , correlated image ellipticities and account for a substantial amount of the total cosmic shear 
^ signal on small angular scales in ACDM and open CDM models. 

^ . 

m ■ 

2 ■ !• Introduction 

. It is now generally agreed that it will be possible in the near future to obtain precision 
' cosmological results via weak lensing measurements. However, this statement is often 
O ^ interpreted to mean merely that with appropriate due diligence on the observational 
' end, precision constraints will result. The focus of this paper is to remind the reader 
^ that highly accurate theory is equally important to the goal of placing strict constraints 
' on cosmology. Simon White expressed this very nicely in his review talk on numerical 
d simulations when he said simply "Precision cosmology will require precision simulations." 
^ ^ This paper attempts to demonstrate the truth of Simon's statement by studying the 
• . effects of multiple weak lensing deflections on both the galaxy-galaxy lensing signal and 
' the cosmic shear signal. Galaxy-galaxy lensing has come a long way since the 1995 
^ . lAU Symposium on gravitational lensing, where it appeared on the programme in the 
■ " " ' "Emerging Techniques" session (Brainerd, Blandford & Small 1996a). Due to a great 
deal of work by a number of different groups, galaxy-galaxy lensing has not only been 
detected with impressively high statistical significance, but it is also being used to place 
strong constraints on the nature of dark matter halos and on the bias between light and 
mass in the universe (see, e.g., Fischer et al. 2000; McKay et al. 2001; Hoekstra, Yee, & 
Gladders 2004; Kleinheinrich et al. 2003; Kleinheinrich 2004; Natarajan et al. 220; Seljak 
2004) 

Simulations of galaxy-galaxy lensing for a moderately deep imaging survey (/um ~ 23) 
performed by Brainerd, Blandford & Small (1996b), hereafter BBS, showed that most 
of the galaxies with magnitudes in the range 22 / ^ 23 would, in fact, have been 
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lensed at a comparable level by two or more foreground galaxies (e.g., §3.5 of BBS). 
That is, multiple deflections are expected to occur quite frequently in deep data sets 
and, in particular, when comparing predictions for cosmic shear from a simulation to 
observations of cosmic shear on small angular scales, it is important that the simulation 
reproduce faithfully all of the weak galaxy lenses along the line of sight. In other words, 
a fair comparison between observations and theory on the scales for which galaxy-galaxy 
lensing is important depends crucially on the ability of simulations to follow the growth 
of the non-linear power spectrum accurately. Even with current codes and computer 
architecture, this is still somewhat challenging and appropriate care must be taken (see, 
e.g., Simon White's review in these proceedings). 

The Hubble Deep Field (North) and the flanking flelds have been the subject of a 
deep redshift survey (Cohen et al. 2000) as well as an extensive multicolor photometric 
investigation (Hogg et al. 2000). As a result, both the redshifts, z, and the rest-frame 
blue luminosities. Lb, of ^ 600 galaxies in this region of space are known (Cohen 2002). 
Therefore, it is possible to make quite a detailed theoretical prediction for the weak 
galaxy-galaxy lensing shear field in the region of the HDF-North and, specifically, for 
the probability and efl'ects of multiple deflections due to galaxy-galaxy lensing. 



2. Halo Properties and Fiducial Model 

For simplicity, the approach of BBS in their galaxy-galaxy lensing simulations was 
adopted and the physical properties of the dark matter halos around the galaxies in the 
HDF-North and the flanking fields were scaled in terms of the characteristic properties 
associated with the halos of L* galaxies. The velocity dispersion of an L* galaxy halo is 
given by cr* and it was assumed that a Tully-Fisher or Faber-Jackson type of relation 
held for each of the galaxies. Therefore, 

4) • 

where is the velocity dispersion of a halo in which a galaxy with luminosity Lb resides. 
The density profile of the galaxy halos was taken to be 

where G is Newton's constant and s is a characteristic halo radius. Further, it was as- 
sumed that the mass- to light ratio of a galaxy was constant independent of its luminosity 
and, therefore, the radii of the halos of galaxies with Lb ^ L*g scale with the radii of the 
halos of L*g galaxies according to: 

The total mass of the halo of an L* lens galaxy is finite and given by: 

M* = (2.4) 

Having made these assumptions, it is then possible to make predictions for the galaxy- 
galaxy lensing shear field within the region of the HDF-North that would be generated 
by the galaxies in Cohen (2002). The lens galaxies include not only galax;ies in the 
HDF-North, but also galaxies in the flanking fields. Because some of the galaxies in the 
flanking fields have quite substantial masses, the weak lensing effects of these galaxies 
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Figure 1. Theoretical shear field in the region of the HDF North. See text for model details. 

have the potential to affect the shear field inside the much smaller region of the HDF- 
North. Therefore, the flanking field galaxies were included in all of the calculations for 
the weak shear field inside the HDF-North itself. Also, in order to have a consistent 
limiting magnitude for the galaxy lenses, only those galaxies in Cohen et al. (2000) with 
R ^ 23 were used in the calculations (i.e., the completeness limit of the full redshift 
survey is deeper in the HDF-North than it is in the fianking fields). 

Shown in Fig. 1 is the theoretical shear field in the region of the HDF-North that 
would be produced by the 427 galaxies in Cohen et al. (2000) and Cohen (2001) for 
which R < 23, and both spectroscopic redshifts and rest frame blue luminosities are 
known. The median redshift of the lenses is ~ 0.6. Source galaxies were assumed to have 
apparent magnitudes in the range 19 < / < 25, and were distributed randomly on the 
sky. A fiat, A-dominated cosmography with Hq = 70 km/s/Mpc, = 0.3, and Aq = 0.7 
was adopted for Fig. 1. A fiducial L* galaxy model for which the velocity dispersion of 
the halo is a* = 150 km/s and the outer scale radius is s* = 100 kpc was used and 
the source galaxies were assumed to follow a redshift distribution of the form: 



P{z\I) 



T{3l(3)zl 



(2.5) 



(e.g., Baugh & Efstathiou 1993), which is in good agreement with the redshift surveys 
of LeFevre et al. (1996) and LeFevre et al. (2004). Here 



Zo = k^[Zm + Z'JJ. - /„)], 



(2.6) 



where Zm is the median redshift, is the median 7-band magnitude, and z'^ is the 
derivative of the median redshift with respect to I. Extrapolating the results of LeFevre 
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et al. (2004) to a sample of galaxies with 19 < / < 25, Zm = 0.86, z'^ = 0.15, = 0.8, 
and /3 = 1.5. The median redshift of the sources is ~ 1.2. Fig. 1 shows the mean over 
6500 Monte Carlo realizations of the weak Icnsing shear field and a close comparison of 
this figure with an image of the HDF-North shows that the peaks in theoretical shear 
field correspond to the brightest galaxies in the HDF-North. 

3. Multiple Deflections 

This probability that a given source will have been weakly-lensed by one or more 
foreground galaxies is, of course, a strong function of the actual value of the shear, 7, 
due to a given weak lensing deflection. That is, it is much more likely for a distant 
galaxy to be lensed by a foreground galaxy which produces an insignificant weak shear 
of 7 ~ 10~^ than, say, a large weak shear of 7 ~ 0.01. Therefore, in order to discuss 
the total number of weak deflections that a given source galaxy is likely to encounter, 
a decision has to be made as to what minimum value of 7 qualifles as a "signiflcant" 
deflection. 

A typical value for the net shear due to galaxy-galaxy lensing is 7 ~ 0.005 (see, e.g., 
the observational papers cited in the Introduction) and this value of 7 was used as a 
baseline for computing the number of weak lensing deflections that source galaxies had 
undergone. Speciflcally, the probability that a given galaxy was lensed by Nd foreground 
galaxies, P{Nu), where each individual deflection gave rise to a shear of 7 > 0.005 was 
computed. That is, P{No — 2) is the probability that a given galaxy had been lensed by 
two individual foreground galaxies, each of which lensed the distant galaxy at a level that 
is comparable to or greater than the expected net shear due to galaxy galaxy lensing. 
Therefore, the results shown here arc an extremely conservative estimate of the frequency 
of multiple deflections. 

Shown in Fig. 2 is the theoretical probability distribution function, P{Ni)), for source 
galaxies in the HDF North with a given redshift, Zg. The fiducial halo model and redshift 
distribution adopted for Fig. 1 were also used here. In addition, the cosmography was 
varied from the flat, A dominated model used in Fig. 1 to include both an open (ilo = 0.3, 
Ao = 0) and an Einstcin-de Sitter model for Fig. 2. From this figure, then, the probability 
that a source at redshift Zs ~ 1 has been lensed at a significant level by two or more 
foreground galaxies is of order 60%. Also, it is clear from Fig. 2 that the probability 
of multiple deflections occurring is essentially independent of the cosmography. Instead, 
the frequency is influenced most strongly by the number of massive lens galaxies that 
are close to the line of sight (i.e., galaxy-galaxy lensing provides information primarily 
about the potentials of the lens galaxies, not cosmology per se; see also BBS). 

At a source redshifts Zg ~ 1.5, the probability of a source galaxy encountering multiple 
deflections of 7 > 0.005 increases to of order 90%. Therefore, in a deep data set for which 
the median redshift is ^ 1, it should be expected that multiple weak deflections of a 
substantial magnitude are very likely to have occurred. Of course, in the case of individ- 
ual deflections for which 7 < 0.005, the probability of comparable multiple deflections 
occurring at any given source redshift will be greater than the results shown in Fig. 2. 

4. Effects of Halo Parameters on Multiple Deflections 

While; the; occurrence of multiple deflections is largely independent of the cosmology, 
it is certainly not independent of the details of the halo parameters. In Figs. 1 and 2 a 
fiducial L* galaxy halo with a* = 150 km/s and s* = 100 kpc was adopted. However, 
galaxy-galaxy lensing constraints on the characteristic velocity dispersion of the lens 
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Figure 2. Probability of a given source galaxy in the HDF-North undergoing Nd weak deflec- 
tions of individual magnitude 7 > 0.005. Panels correspond to sources with different redshifts, 
Zs- Solid histogram: flat, A-dominatod model. Dashed histogram: open model. Dotted histogram: 
Einstein-de Sitter model. Vertical dotted line shows No = 2. 



galaxies range from a* ^ 135 km/s (e.g., McKay et al. 2001; Hockstra, Yee, & Gladdcrs 
2004) to (7* ~ 165 km/s (e.g., Kleinheinricli ct al. 2004). Constraints on the characteristic 
scale radius are few, and based on the galaxy galaxy tensing observations of Hoekstra, 
Yee, & Gladders (2004), the halos of L* galaxies may be as large as s* ~ 185 kpc. 

In this section, then, the effects of varying the halo parameters on the occurrence of 
multiple deflections is investigated. Since the occurrence of multiple deflections is only 
weakly dependent on the cosmography, for the remainder of this paper only the flat, A- 
dominated model was used. Again, all source galaxies were assumed to have magnitudes 
in the range 19 < 7 < 25, with redshifts determined by eqns. (2.4) and (2.5) above. 
The halo parameters of the lens galaxies in the HDF-North and flanking fields were then 
varied as follows: a* = 135 km/s, 150 km/s, and 165 km/s; s* = 50 kpc, 100 kpc, 
and 200 h''^ kpc. 

Shown in Fig. 3 is P{Nd) for all source galaxies in the HDF-North with 19 < 7 < 23. 
Like Fig. 2, the minimum value for a "deflection" to be counted in this figure is 7 = 0.005. 
Unlike Fig. 2, however, P{Nd) has been computed over the entire redshift distribution of 
the sources. From this figure, then, the probability that a given source in the HDF -North 
has been tensed more than once is 27% for the lowest mass L* halo (upper left panel), 
53% for the fiducial halo of Figs. 1 and 2 (center panel), and 69% for the highest mass 
L* halo (bottom right panel). Therefore, it is clear that the mass adopted for the halo of 
an L* galaxy has a rather substantial effect on the number of multiple deflections that 
occur in a galaxy-galaxy lensing calculation. 

Since the halos are roughly isothermal and, therefore, the shear decreases with pro- 
jected radius approximately as ^{6) oc 0~^, it is interesting to ask whether the closest 
tens to a given source (in projection on the sky) is necessarily the strongest tens. That is, 
for a source that undergoes multiple deflections, is the strongest deflection most likely to 
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Figure 3. Probability of a given source galaxy in the HDF Nortii undergoing No weak deflec- 
tions of individual magnitude 7 > 0.005. Panels correspond to different halo models, indicated 
by the values of a* and s* . In all cases a ACDM model was used. Unlike Fig. 2, here P{Nd) 
was computed over the entire redshift distribution of sources with 19 < / < 25. Vertical dotted 
line shows No = 2. 



come from the nearest lens on the sky? The answer to this question is "No," and is shown 
clearly by Fig. 4. In the case of sources for which the net shear (i.e., the total shear after 
all weak deflections have occurred) is 7 ^ 0.01, more than 50% of the time the strongest 
lens is not the closest lens on the sky. Note that this result is essentially independent 
of the characteristic halo parameters that were chosen, which is consistent with the fact 
that, by and large, it is only the most massive galaxies that contribute substantiafly to 
the overall galaxy-galaxy tensing signal. 

Figs. 5 and 6 demonstrate a result that is perhaps somewhat counter-intuitive: multiple 
weak tensing deflections, on average, give rise to a greater net shear on the source galaxies 
and a greater mean tangential shear about the lens centers. That is, multiple deflections 
do not simply cancel one another out, leading to little or no net shear on the source 
galaxies. Fig. 5 shows probabilities for the distribution of the ratio of the maximum 
value of any given single deflection, 7max, to that of the net shear, 7net, for alt sources 
with 19 < 7 < 25 that underwent more than one deflection of any magnitude (i.e., 
7 > 0). The vertical dotted tine shows the median value of the distribution, and for all 
halo models the median is less than 1. That is, more than 50% of the time, the net shear 
on a given source is greater than the maximum single deflection that it underwent. 

Fig. 6 shows the mean tangential shear, computed about the lens centers in the HDF- 
North. Squares show the results of the proper inclusion of multiple deflections for all 
sources, and crosses show the result of tensing each source solely by the closest tens on 
the sky (i.e., a "single deflection" calculation). In all cases, the inclusion of multiple 
deflections gives rise to a larger mean tangential shear on scales 6 ^ 1", and in the case 
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of the more massive halos, the increase in the shear is quite substantial. This figure, then 
emphasizes the need for a correct, multiple deflection calculation when using observations 
of the mean tangential shear to constrain the halo properties of the lens galaxies (i.e., 
the comparison of single deflection calculations to observations leads to an inferred halo 
mass for L* galaxies that is too large). 



5. Correlated Image EUipticities and Cosmic Shear 

In addition to giving rise to a generally larger net shear on source galaxies and a 
larger mean tangential shear about the lens centers, multiple deflections in galaxy-galaxy 
lensing give rise to correlated ellipticities in the images of the galaxies. This is, of course, 
precisely the eflfect of cosmic shear, but in the case of galaxy-galaxy lensing, this is merely 
the very small k end of the power spectrum (i.e., the highly non-linear regime) that is 
contributing to the total cosmic shear signal (i.e., when computed over all structures 
along the line of sight). 

The degree to which multiple deflections give rise to correlated image ellipticities and, 
hence, contribute to the cosmic shear signal is a strong function of the mass adopted for 
the halo of an L* galaxy. This is shown in Fig. 7, where the shear correlation function, 

C^^(0) s (7i • 7*) , i^j (5.1) 

is shown as a function of the halo parameters. The mean value is computed for all 
foreground-background pairs of galaxies separated by angles 9 ± 66/2 on the sky (see, 
e.g., Blandford et al. 1991). Here 7i is the image shape of galaxy i and jj is the complex 
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Figure 5. Probability distribution for the ratio of the maximum shear due to a single deflection, 
7max, to that of the the net shear, 7nct- Vertical dotted line shows tlic median of the distribution. 
All sources with 19 < 7 < 25 that underwent more than one deflection of magnitude 7 > were 
included in the calculations. 

conjugate of the image shape of galaxy j. The image shape is defined as 



7; 



a2 + 62 



(5.2) 



whore a and b are the major and minor axes of the image equivalent ellipse and (j) is 
its position angle. As in Fig. 6, squares show the results of the full multiple deflection 
calculations and crosses show the results of calculations in which each source was lensed 
solely by the closest lens on the sky. From this figure, then, the single deflection cal- 
culations do not give rise to correlated image ellipticities on scales 9 ^ 5" (on smaller 
scales, of course, the images of distant sources that have been lensed by the identical 
foreground galaxy will be correlated because of the tangential alignment about the lens 
center). In the multiple deflection calculations, however, sufficiently massive halos give 
rise to correlated ellipticities that persist to significantly large angles, just due to the 
galaxy-galaxy lensing signal alone. 

Oftentimes the phrase "cosmic shear" is interpreted to mean "lensing by large-scale 
structure" , and this is true on large angular scales where it is only structure in the linear 
regime that is responsible for the gravitational lensing. Properly, however, cosmic shear 
is the lensing of distant galaxies by all the mass along the line of sight, including highly 
non-linear structures. This is why, on small angular scales, it is necessary to use large 
simulations (e.g., Jain, Seljak, & White 2000; Valageas, Barber, & Munshi 2004; Vale & 
White 2003) to make accurate theoretical predictions of cosmic shear. 

Prom Fig. 7, it is clear that, depending upon the characteristic mass of the halos of L* 
galaxies, galaxy-galaxy lensing will contribute to the cosmic shear signal as measured, 
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Figure 6. Mean tangential shear computed about the lenses in the HDF-Nortli. Squares show 
the results of the multiple deflection calculations, crosses show the result of lensing each source 
gala^xy solely by the nearest lens on the sky. All sources with 19 < 7 < 25 were included in the 
calculations. 



for example, via the top hat shear variance, 

2 f°° dk 



p^ik) [Jiike)Y 



(5.3) 



Here Pk is the power spectrum of the projected mass density of the universe, Ji is a 
Bessel function of the first kind, and 6 is the size of the aperture over which the mean is 
computed. In an observational data set, the function is computed as 



N{N-1) 



(5.4) 



for all galaxies within an aperture of size 9 on the sky (see, e.g., Hammerle et al. 2002). 
The above published predictions for (7^) in CDM models have been based on simulations 
in which the mass and force resolution were not quite adequate to resolve the halos of L* 
galaxies particularly well and, so, it is difficult to be certain how accurate the theoretical 
predictions truly are on very small angular scales. 

1/2 

The symbols in Fig. 8 show the r.m.s. cosmic shear, (7^) , computed within circular 
apc;rtiu-es of radius 6 in the full multiple deflection calculations of galaxy-galaxy lensing 
in the HDF-North. That is, the symbols in this figure show the r.m.s. cosmic shear due 
to galaxies alone, without any contribution from linear or quasi-linear structures along 
the line of sight. Ftom bottom to top, the point types correspond to increasing the halo 
mass from a minimum mass of M* — 0.83 x lO^^M© (open diamonds) to a maximum 
mass of M* = 4.05 x lO^^M© (four-pointed stars). The mass of the fiducial halo with 
a* = 150 km/s and s* = 100 h'^ kpc is M* = 1.67 x lO^^M©, and the r.m.s. cosmic 
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Figure 7. Shear correlation funetion, C^f-,{9), as a funetion of angular seale and halo parameters. 
Squares show the results of the multiple deflection calculations, crosses show the result of lensing 
each source gala^xy solely by the nearest lens on the sky. All sources with 19 < / < 25 were 
included in the calculations. 



shear in this case is shown by the solid squares. Also shown is the prediction for (7^) 
for ACDM and open CDM from the simulations of Jain, Scljak, & White (2000) for a 
source galaxy population with median redshift of ^ 1.2 (i.e., similar to the HDF- 
North raytracing simulations for sources with 19 < / < 25). The ACDM and OCDM 
predictions do not extend below 9 0.4'. 

The filled circles and filled triangles show the r.m.s. cosmic shear produced by galaxy- 
galaxy lensing alone for cases in which the mass of the halo of an L* galaxy is ~ 20% less 
than the fiducial halo mass and ~ 20% greater than the fiducial halo mass, respectively. 
At 6* ~ 0.5', then, a change in the fiducial halo mass of only 20% results in a change 
in the predicted cosmic shear signal that is quite comparable to the predicted differ- 
ences between two rather different cosmological models. In other words, if one wishes to 
use observations of cosmic shear on small angular scales to constrain the cosmography 
(through, e.g., a comparison to ray-tracing simulations) it is vital that the simulations 
have followed the formation of the growth of non-linear structures very accurately. 

The r.m.s. cosmic shear due solely to galaxy-galaxy lensing extrapolates to zero at 
^ ~ 1' for the fiducial halo model with a* = 150 km/s and s* = 100 kpc. It is, 
therefore, only on scales <, 1' that contributions of galaxy-galaxy lensing to the cosmic 
shear signal are likely to be of importance. 



6. Conclusions 

The occurrence and effects of multiple weak deflections due to galaxy-galaxy lensing 
were investigated for a deep data set in which 2;iens ^ 0.6 and -Zsource ^ 1-2. Ray-tracing 
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Figure 8. The r.m.s. cosmic siiear as a function of aperture size, 8. Points show the results 
from ray-tracing simulations of galaxy-galaxy lensing in the HDF-North for which all sources 
with 19 < / < 25 were included in the calculations. Prom bottom to top, the points correspond 
to increasing halo mass (see text). Filled squares show results for a fiducial L* galaxy halo for 
which a* = 150 km/s and s* = 100 kpc. Filled circles and filled triangles show results for 
a ~ 20% decrease in the fiducial halo mass and a ^ 20% increase in the fiducial halo mass, 
respectively. Dashed and dotted lines show theoretical predictions for open CDM and ACDM, 
respectively, based on calculations from Jain, Seljak, & White (2000). 



simulations of galaxy-galaxy lensing by lens galaxies with i? ^ 23 in the Hubble Deep 
Field (North) and flanking fields were used to compute the net shear on source galaxies 
with magnitudes in the range 19 < / < 25. Both the redshifts and the rest frame blue 
luminosities of the lenses are known, which allows for a detailed theoretical prediction of 
galaxy-galaxy lensing in the HDF North, given a particular cosmography and a model 
for the halos of the lens galaxies. The primary conclusions from this work are: 

• Multiple weak deflections are commonplace in such deep data sets, and the proper 
inclusion of multiple deflections is important to a correct prediction for the net shear 
experienced by the majority of source galaxies. 

• For a given source redshift, the probability of multiple weak deflections is largely 

insensitive to the cosmography (i.e., galaxy galaxy lensing is miich more sensitive to the 
details of the gravitational potentials of the lens galaxy halos than it is to the values of 
the cosmological parameters). 

• Compared to a single deflection calculation in which sources are lensed solely by the 

nearest lens on the sky, a full multiple deflection calculation leads both to a larger net 
shear for most individual sources and a larger tangential shear about the lens centers. 
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• Multiple weak deflections give rise to correlated image ellipticities (i.e., "cosmic 
shear" due to power on highly non-linear scales). 

• On angular scales 9 < I' , galaxy-galaxy lensing alone accounts for a substantial 
amount of the cosmic shear signal expected in ACDM and open CDM models. The 
magnitude of the signal is, however, very sensitive to the characteristic halo mass and, 
hence, accurate comparisons of observations of cosmic shear and theory on such scales 
rely heavily on the ability of simulations to follow the growth of the non-linear power 
spectrum with high accuracy. 
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